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Abstract —Fifth generation networks will incorporate a variety 
of new features in wireless networks such as data offloading, 
D2D communication, and Massive MIMO. Massive MIMO is 
specially appealing since it achieves huge gains while enabling 
simple processing like MRC receivers. It suffers, though, from a 
major shortcoming refereed to as pilot contamination. In this 
paper we propose a frame-work in which, a D2D underlaid 
Massive MIMO system is implemented and we will prove that 
this scheme can reduce the pilot contamination problem while 
enabling an optimization of the system spectral efficiency. The 
D2D communication will help maintain the network coverage 
while allowing a better channel estimation to be performed. 

I. Introduction 

The rapid increase in smart-phones number allowed the 
emergence of new services that although, enriching for the 
user experience, lead to an exponential growth in the traffic 
volume over the wireless networks. Since the network re¬ 
sources especially the used spectrum, are limited, techniques 
like Massive MIMO, D2D communication, data offloading and 
small cell networks were proposed in order to provide the 
necessary increase in network capacity for future generation 
networks. In addition to traffic volume increase, the issue of 
energy efficiency of the network has to be dealt with. To this 
end two of the most promising techniques are Massive MIMO 
and D2D underlaid network. 

Multiuser Massive MIMO systems were intensively in¬ 
vestigated and it was proven that with simple processing 
this systems allow an increase in network capacity with the 
possibility of a more energy efficient network. Nevertheless, 
Massive MIMO systems are no without shortcomings. In fact 
while the large number of antennas mitigate the effect of 
uncorrelated noise and fast fading 0, the interference caused 
by pilot reuse persistf2|. The pilot contamination problem was 
addressed several times, in 0 a subspace projection method 
was used to mitigate the pilot contamination problem.In 0, 
the authors proposed a covariance-aided channel estimation 
framework and a coordinated pilot assignment strategy in order 
to reduce pilot contamination impact. Another promising tech¬ 
nology is Device to Device communication. This technology 
allows the user equipment to communicate directly using the 
cellular network resources instead of going through the cellular 
infrastructure. 

With enabled D2D communication, proximity-based ser¬ 
vices become possible 0. This proximity enable an increase 


in spectral and energy efficiency. It make also more efficient 
data offloading schemes possible. Combining D2D commu¬ 
nication and Massive MIMO, has been proposed in m and 
IH for MIMO networks in FDD mode and has proven to 
enable more efficient forms of feedback. In this paper, we 
investigate the coexistence of the two technologies in TDD 
mode. In order to reduce the pilot contamination phenomenon, 
we consider in addition an exclusion zone framework and 
develop an optimization framework that determines the best 
exclusion zone size and power allocation. 

II. Contribution and outcomes 

The contributions of our work are summarized as follows: 

1) A base station controlled mode selection scheme: 

In this paper we propose a mode selection scheme based 
on the user’s position. This scheme divide the cell into 
two parts where cellular mode and D2D mode are not 
allowed to coexist|14j- The Base station will fix the 
radius of the exclusion zone in a scalable way based on 
the user density, user transmit power, average cellular 
SINR , and average interference on D2D links. 

2) Impact of the mode selection on Pilot contamination: 

In this paper we consider a realistic framework in 
which the Base station uses Pilot-based CSI estimation 
only for cellular users. This enable us to assess the 
impact of our proposed mode selection scheme on a 
major shortcoming of Massive MIMO systems which 
is Pilot contamination. We And that the mean square 
error of channel estimation due to Pilot reuse is reduced 
in our system but vary depending on whether D2D 
communication is allowed or not during training phase. 

3) Average SINR Optimization: 

We develop the analytical expressions of the average 
SINR for a reference cellular user and the average inter¬ 
ference on a reference D2D link. Using this expressions 
we introduce a quasi-concave optimization problem in 
which cellular user power and the exclusion region 
radius are defined so that the average cellular SINR is 
maximized while an acceptable interference level on the 
D2D communication is guaranteed. By using the average 
interference and average SINR we are able to perform 
optimization without the need for full CSI of all users. 


This is quite practical in a dense system like the one 
considered in this paper. 

III. SYSTEM MODEL AND PRELIMINARIES 

In this paper we consider a Multi-user multi cell MIMO 
system with a large excess of base station antennas operating 
in time-division duplex (TDD) mode in which the coherence 
time of the channel is divided between a training and a 
data transmission phase. The base station will use channel 
reciprocity and a mutually orthogonal set of pilots in order to 
estimate the channel of each scheduled cellular user. 

We consider also that in addition to the cellular users the 
network is underlaid with D2D terminals using the same 
frequency spectrum as the cellular users. 



Fig. 1: Network model 

The cellular and D2D users locations are modeled as a 
Poisson Point Process (PPP) denoted <I> with density A . In our 
framework the users will select either cellular mode or D2D 
mode based on their distance from the serving Base station. 
The Base stations follow also a PPP <l>/,with density A/,. Each 
Base station will be serving users within a cell radius of R c 
and each Base station will be equipped with M antennas. 
The D2D and cellular users transmit at powers P,i and P c 
respectively. We define P c = CPd where C is a parameter 
calculated in order to optimize the system performance and 
A = a x At,. Note that in this paper we always consider 
P c > Pd . In fact.while D2D links are more suited for short 
range communication the cellular links are utilized for longer 
range. In our setting priority is given for the cellular users since 
they are also expected to play the role of a relay point for D2D 
user wishing to communicate through the infrastructure of the 
cellular network, so the quality of the link of cellular user 
needs to be high. We can use, instead, some low power access 
points at the edge of the exclusion zone in order to relay 
data coming from D2D users to the cellular infrastructure. 
This setting enables us to guarantee the accessibility to the 
cellular infrastructure using low power links which should 
increase the energy efficiency of the network and reduce the 
interference. We will focus on the central Base station located 
at the origin. The received signal at this base station during 
data transmission phase will be: 
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Where a > 2 denotes the pathloss exponent, r k b and r, 
denote, respectively, the distance of the kth cellular user in 
the bth cell and the ith D2D user from the reference Base 
station. h kb £ (jMxi anc | e (jMx i respectively, 

the vector channels from the kth cellular user in the bth cell 
and the ith D2D user to the reference Base station with i.i.d 
CN(0 , l)elements. u k b and tq denote the zero mean unit 
variance transmitted data symbols of the kth cellular user in 
the bth cell and the ith D2D user. N c is the additive white 
Gaussian noise at the reference Base station. 

The received signal at a reference D2D receiver in the 
reference Base station during the uplink phase is given by: 
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Where rjkb and r :jl denote, respectively, the distance of the 
kth cellular user in the bth cell and the ith D2D user from the 
jth D2D receiver taken here as a reference. g ;) kh eC lxl and 
gji £ C lxl denote, respectively, the vector channels from the 
kth cellular user in the bth cell and the ith D2D user to the 
jth D2D receiver. N d £ C lxl is the additive white Gaussian 
noise at the jth D2D receiver. Gaussian signaling is assumed 
and all channel coefficients are i.i.d C'A'(0.1) . In this paper 
we consider MRC receivers which proved to be sufficient to 
achieve Massive MIMO gain. We will focus on the imperfect 
CSI case where only cellular users channel is estimated. The 
idea developed here is to restraint the location of the cellular 
users within the cell which will reduce on average the impact 
of pilot contamination. In fact, while in the classical network 
co-pilot users may be located anywhere within the cell, in our 
framework these users will be within a circle of radius Re 
centered on the base station and thus increasing the average 
distance between co-pilot users. We consider the Matrix Q £ 
C TpxT p containing the pilot sequences that satisfies Q*Q = 
It the same pilot sequences are reused in every cell with a 
reuse factor of 1. In this case the M xT p dimensional received 
signal at the reference Base station during the training phase 
is: 

y c = EE \fPc\\Tkb\\ 2 hkbqk* + 
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with vF i represent the T p x 1 data vector sent by the 
D2D transmitters during the training phase, NJ! the M x T p 
dimensional additive Gaussian noise at the reference Base 


station and q k £ C TpXl the pilot vector used by user k in 
cell b. 

IV. Average cellular SINR & average 
INTERFERENCE ON D2D RECEIVERS 

Upon receiving the training signal, the Base station will 
estimate the channel of the kth cellular user in its cell as 
follow: 

hki = YcQk 

The channel estimate in this case depends on whether D2D 
communication is allowed or not during training phase. In 
fact while allowing D2D transmitters to be active during this 
phase may allow us to use the system resources efficiently it 
may on the other hand introduce a further degradation of the 
channel estimate since, in addition to the usual limiting factor 
of pilot contamination, further degradation due to active D2D 
transmitters persist even in a Massive MIMO systemfT]. 

Two expressions of the estimate of the channel of kth user can 
be derived: With muted D2D the channel estimate will be: 

hki = ^ 2 VP~c\\r k b\\~ h kb + NPq k 

b 

If D2D communication was active during training phase the 
channel estimate is: 

h k i = ^2 VK\\rkb\\^ hkb+^2 VPdWnW^ hiU^*q k +NPq k 

b i 

A. Interference analysis 

In the framework of our system the reference base station 
will be subject to interference coming both from D2D and 
cellular users. The cellular interference will be generated by 
co-pilot users in other cells. This is the usual limiting factor 
of Massive MIMO system the rest of the cellular interference 
vanish thanks to the large number of receive antennas at the 
Base station. Since the cellular interference results from the 
pilot contamination phenomenon, with a pilot reuse factor of 
1 the cellular interferer will have the same density X b as the 
Base stations. 


radius D. Then, the remaining points of (f> n form a Poisson 
hole process with density \ n = X n exp (— A m 7rD 2 ). 

The new density of the D2D users is then given by: 

A d = A exp (—7 tA b Re 2 ) 

for a D2D receiver the interference will come from all 
cellular users and from the other D2D transmissions.We con¬ 
sider that any D2D receiver only knows the channel to its 
transmitter. 

for a typical D2D receiver during uplink the cellular in¬ 
terference will come from cellular user spanning over circles 
with radius R e centered on the network Base stations. Then 
for the D2D receiver the cellular interferer will have a density 
of: 

A c = A(1 - exp (-7 tA b R 2 e )) 

while the density of the D2D interferer is A d 

B. Impact of the exclusion zone on channel estimation 

In order to show the impact of the exclusion zone on the 
CSI acquisition we calculate the mean square error of channel 
estimation. We can show then, that with the exclusion zone 
and deactivated D2D users during the training phase, we are 
able to reduce the Mean square error in channel estimation 
we also can see that this error increases when increasing the 
radius of the exclusion zone, the mean square error of the 
channel estimation is given by: 

MSE = E(e c hannel£ c hannel) = E((h k i — h k \)(h k \ — h k i) ) 

with deactivated D2D transmission during training phase the 
MSE after normalization is: 

MSE = BUY, IMI^ h kb + Nfq k - h kl f) 

b 
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If the D2D transmission is active during the training phase 
the MSE will have an additional term coming from the D2D 
interference which will reduce the accuracy of the channel 
estimation. In this case the MSE becomes: 


D2D interference will persist if D2D transmission was 
active during the training phase, but the interfering signal 
will come from limited areas defined by the exclusion area 
radius. In fact in each cell the D2D terminals are allowed 
to be active only within a ring of inner radius Re and outer 
radius R c . Limiting the area on which D2D interferer may 
span changes the density of active D2D users since under our 
system assumptions the D2D users form now a point process 
referred to as Poisson Hole Process which is a Cox process 
0 . 

Definition 1 (Poisson hole process): 

let 4> m and rj> n two PPPs with respective intensities A m and 
X n with A„ > X m . For each x £ fm, remove all points in 
(f> n fl b(x, D) where b{x , D) is a ball centered at x and with 


MSE = Mi — 
y Pc 


2irXb{2R c - R e ) 2 ~ a 
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P d 2nX d {R e ) 2 ~ a 
(a - 2 )P C ’ 


C. Average cellular and D2D SINR 

In order assess the performance gain generated by the 
exclusion region framework we will investigate the average 
SINR of a reference user located at the reference cell. With 
the assumption that this user (k) is in cellular mode we have 
after receiving the uplink signal and using a MRC receiverfTI: 
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Conditioned on the position of user k taken here as a 
reference user we can have the following average SINR 
expression: 


E(SINR kl \r kl ) 
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Proof: See Appendix A. 


We notice the persistence of the D2D interference even with 
an infinite number of antennas. This due to the fact that D2D 
interference is active during the training phase m . One way to 
remove the D2D interference is to mute D2D communication 
during training phase.In this case the average reference user 
SINR will become: 


cvoriUD L, t ll 7 ’ fc l|| 
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For computing the average SINR of a typical D2D receiver 
we consider a user in D2D mode located at a distance d 
from the reference base station and having a distance r :j from 
its transmitter, for a typical D2D receiver the average SINR 
conditioned on its position is : 

E(SINR d j \r j ,d ) 

= _ P d \\d\\~ a _ 

_2 I P c 27rA c (d-Pe) 2 -“ , P d 2-K\ d (r 0 ) 2 -° L 
a D2D ~r a —2 a-2 

V. Optimizing the cellular user average SINR 

Since the major shortcoming of D2D underlaid cellular net¬ 
work is the cross tier interference. Optimizing the respective 
cellular and D2D transmission powers and the exclusion zone 
radius should enable us to extract the full potential of our 
system setting. We propose an optimization problem in which 
transmission power and the radius of the exclusion zone are 
optimized. We concentrated on the average interference and 
average SINR because we don’t need perfect CSI in this case. 
In fact while the usual Power Control problems require perfect 
knowledge of the user’s channels, our problem only deals 
with the average interference which can be computed without 
exact knowledge of the channels. This is quite practical, since 
acquiring channel information of all active users, both cellular 
and D2D, in the cell will cause a reduction in the data 
transmission phase. 

Only basic parameters of the network are needed to find the 
optimal solution in our problem. Taking X = ( R ei C ) as the 
vector of optimization variables, f(X) = E(SINR re f\r re f) 
and g(X) = P^N(d-Re) 2 - + p d 2.x d (r 0 ) 2 - we can formu _ 

late the optimization problem as: 

maximize f(X) 
subject to g(X) < Id 2 D 

Our optimization problem enable us to couple the exclusion 
region radius optimization with a power control problem 


allowing effective interference mitigation both for the D2D 
and cellular communication. Although the optimization 
problem is not convex, the solution can be derived with zero 
duality gap. 


lemmal 1. The optimal solution of the proposed optimization 
problem is found using KKT CONDITIONS. And it is given 
by the two following equations: 

4Cn\ d \\r kl \r 2a (R e ) 2 ~ 2a on 2nX c (d - Re) 2 ~ a _ n 
D 2 (2a — 2) 1 d a-2 

- ||rfct|r 2 “ [(2a - 2)\ b C\2R c - R e ) l ~ 2a + A 
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x (C((a 2)(d R e y~ a A c + 2§ 

(d - R e ) 2 ~ a A exp (-^) 2 ) - A d r 2 -“2^f)] = 0 
where D = g 2 ^A b ( 2 R c -R e ) 2 ~ 2 ° 2 *\ d (R e ) 2 - 2 °‘ md ^ ig f/jg 
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Lagrangian multiplier. 

Proof: See Appendix B. 

■ 

VI. Numerical Results 

In this section we present some numerical results to validate 
the derived analytical model and to investigate the impact of 
the exclusion zone frame-work on the system performances 
depending on the system parameters. First we need to validate 
our analytical model. To this end we compare the simulation 
results for the average received interference at the Base 
station and mean square error in CSI acquisition with the 
results derived from the developed expressions of these two 
quantities. In our simulations we consider 31 hexagonal cells 
with a side length Rc = 1 Km. The resulting Base station 
density can be approximated by A& = gigi for user density 
we consider A = a A/,. The exclusion region radius R e will 
vary from [Remin, Reman]- We suppose that the reference 
user will be located at a distance d = 200m from its serving 
Base Station. 


a 

3 

A b 

I 

7T.R2 

R c 

1 Km 

A 

150A 6 

Re 

[OAKm, 0.9 Km] 

a" 

10"3 


TABLE I: Simulation Parameters 


In the simulation, we generate mobile users according to 
a PPP with density A These users will be randomly located 
in the network. Transmit power will be allocated to users 
according to there locations. The simulation scenario consider 
pilot aided channel estimation. The pilot reuse factor is 1. After 
channel estimation, data detection is performed. The received 
signal at the base station is averaged after 10000 iteration and 
measurements will be done for different values of R e . 
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Fig. 2: Average received interference at the reference BS 


Fig. 4: Average cellular user SINR for different users densities 


We can see that the analytical and simulation results match. 
The two curves represent respectively the analytical average 
received interference at the base station given by the denom¬ 
inator of the cellular user average SINR and the measured 
value after simulation. 

For the MSE we compare the simulation curve with the 
analytical expression and with the MSE curve when no exclu¬ 
sion zone is implemented. When R e = II, i.e the cellular user 
can be anywhere within the cell, the MSE curve represent the 
classical case without exclusion zone. We notice that using our 
system setting the MSE of channel estimation remains below 
the MSE of the classical case and they are only equal when 
R e = R c . 



Fig. 3: Normalized Mean square Error in reference user 
channel estimation 

The next figure shows the average SINR curve evolution as 
a function of R e for different user densities. User density is an 
important parameter in the system. In fact while the cellular 
interferer always have as density A & , the D2D interferer 
density increase while increasing A 

We present, in table II, the optimal values of the opti¬ 
mization variables for different interference thresholds on a 
reference D2D receiver located at 0.95Km from the reference 
Base Station for a D2D transmission power of 16 dBm and a 
Cmax — 10 • 


Id2d(W) 

12 

15 

18 

R e (Km) 

0.7465 

0.5994 

0.6924 

C 

2.2176 

8.6196 

10 


TABLE II: Solutions for different Ijj 2 D 


VII. Conclusion 

In this paper, we investigated the interplay between Massive 
MIMO and D2D underlaid cellular network where the mobile 
users choose their operating mode depending on their location. 
Based on the exclusion zone framework we investigated the 
average user SINR and the mean square error in channel 
estimation. We have shown that this system setting enables 
the reduction of the impact of pilot contamination. We also 
developed an analytical model that allow the optimization of 
the exclusion region radius and cellular user power in order 
to optimize the network performance. Although, the proposed 
framework allow the coexistence of cellular mode and D2D 
mode while minimizing the inter-tier interference, the fairness 
criterion should be investigated especially in terms of energy 
consumption. 


VIII. Appendix 


A. proof for average SINR 

The average D2D interference at the reference Base Station 
is given by: 
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The integral is due to Campbell’s formula 0- In a similar 
way we can derive the expression of the average cellular 
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B. Proof of Lemma 1: 

To prove Lemma 1 we use two propositions given by Arrow 
and Enthoven in Q for quasi-concave optimization problems. 
First we need to prove that the interference to which a typical 
D2D receiver is subject is a quasi convex function of the vector 
of the optimization variables and that the Objective function is 
also quasi concave as a function of the optimization variables. 

In order to prove the strict quasi-concavity of the objec¬ 
tive function it is sufficient to prove the strict concavity of 
'T(X) = Log(f(X)) J8] We compute the hessian Hq ,of 
vP(X) = Log(f(X)) 


Hy = V 2 T'(X) 


taking \ b = —Lj, P c = C x Pd and A = aX b we have: 

7TiX c 

f (X) = Log(C 2 R 2 c (a - 1) ||r fcl ||- 2Q )- 
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In addition we have gg^ > and gg^ > WC 

Then //,[, is a negative semi-definite matrix which proves the 
strict quasi-concavity of f(X) 

For the constraint g(X) it is sufficient to show that the 
interference is an increasing function of the exclusion radius 
R e in fact 


dg 

dR e 

x (C((a - 2 )(d - R e y~ a (l exp (-^) 2 )+ 

a — 2 K c 

2 |( d - ft)2 - 0 exp(-|n- r -2|) 

we can see that this derivative is positive in the interesting 
regime where P c > Pd 

> 0 iff C > C min 


Since X7g(X) > 0 , g(X) is a quasi-convex function 


Now, we use the following propositions: 

Proposition 1 (Arrow and Enthoven): Assume that gi, ■.■,g m 
are quasi-concave functions and that the following 
regularity conditions hold: 

1) Slater condition is verified : 3Xo £ 

C 2 suchthat gi(X 0 ) > 0 V* 

2) for each i g, t is either concave or V<?i(x) Y Ofor each 
feasible solution of the optimization problem 

Then if X* is a local optimal solution of the optimization 
problem then 37 * such that with (X*,y*) the Kuhn-Tucker 
conditions hold. 

Proposition 2 (Arrow and Enthoven): If in addition the 
objective function f is twice differentiable on the feasible set 




























and additionally V/(X*) ^ 0 holds then X* is an optimal 
solution of the optimization problem. 


Cellular Networks, proc. of IEEE PIMRC 2010, Istanbul, 
Turkey, Sep. 2010. 


It is clear that our optimization problem satisfies the opti¬ 
mality conditions stated in the two propositions which means 
that the optimal exclusion zone radius and cellular power can 
be found using KKT conditions. 
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